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ABSTRACT 

The Differential Emission Measure (DEM) analysis is one of the most used diagnostic tools 
for solar and stellar coronae. Being an inverse problem, it has limitations due to the presence of 
random and systematic errors. We present in theses series of papers an analysis of the robustness 
of the inversion in the case of AIA/SDO observations. We completely characterize the DEM 
inversion and its statistical properties, providing all the solutions consistent with the data along 
with their associated probabilities, and a test of the suitability of the assumed DEM model. 

While Paper I focused on isothermal conditions, we now consider multi-thermal plasmas and 
investigate both isothermal and multithermal solutions. We demonstrate how the ambiguity 
between noises and multi-thermality fundamentally limits the temperature resolution of the in- 
version. We show that if the observed plasma is multi-thermal, isothermal solutions tend to 
cluster on a constant temperature whatever the number of passbands or spectral lines. The 
multi-thermal solutions are also found to be biased toward near isothermal solutions around 1 
MK. This is true even if the residuals support the chosen DEM model, possibly leading to erro- 
neous conclusions on the observed plasma. We propose tools to identify and quantify the possible 
degeneracy of solutions, thus helping the interpretation of DEM inversion. 

Subject headings: Sun: corona - Sun: UV radiation 

1. Introduction sure of the amount of emitting plasma along the 

Line Of Sight (LOS) as a function of the electron 
A convenient approach to study the ther- temperature T e . However, the intrinsic difficul- 

mal structure of the solar and stellar outer at- tieg i nvo i vec j m this inverse problem lead to many 

mospheres is the Differential Emission Measure complications in its inference, making its intcr- 

(DEM) formalism. The DEM £(T e ) is a mea- pretation ambiguous. The central point of these 
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series of papers is to provide new tools to sys- 
tematically and completely characterize the DEM 
inversions and assist the DEM interpretation. Us- 
ing the technique develope d for this purpose, ex- 



haustively described in ([Guennou et al. I [2012, 



hereafter Paper I), it is possible to determine and 
to compare the DEM diagnostic capabilities of 
given instruments. With only three bands, the 
previous generation of imaging telescopes were 
shown not to be wel l-suited to DEM analysis (e.g. 
Schmelz et al.|[2009h . The situation has changed 
with the availability of new multi-band instru- 
ments such as the Atmospheric Imagin g Assembly 
(AIA) telescope ( Lemen et al. I l2012t ) on board 
the Solar Dynamics Observatory (SDO). Apply- 
ing our technique to AIA, we show the increased 
robustness of the inversion for isothermal plasmas, 
but we also found intrinsic biases if the observed 
plasma is multithermal. 

Using the notations of Paper I, the DEM is de- 
fined as 



t(T e )=nl(T e )dp/dQogT e 



(1) 



where v? e is the square electron density averaged 
ove r the portions dp of the LOS at temperature 
T e (jCraig fc Brownlll976l ). The observed intensity 
in the spectral band b of an instrument can be 
expressed as a function of the DEM £ as follows 



h = ^~ f R b (n e ,T e )£{T e )d\ogT e 
4tt Jo 



(2) 



Details and references about the DEM formalism 
can be found in Paper I. Given a set of obser- 
vations in N different bands, the DEM can be 
in principle inferred by reversing the image ac- 
quisition process described by Equation [2l How- 
ever, solving for the DEM proved to be a con- 
siderable challenge. The complications involved 
in its derivation are one of the reasons for con- 
troversial results regarding the thermal structure 
of the corona. For example, the DEM is one 
of the method used to derive the still debated 



physi cal properties of plumes (jWilhelm et al 
20111 ). Also, while heating models predict dif- 



ferent thermal structures for coronal loops de- 



pend i ng on the processes i nvolved (e.g. iKlimchuk 



20061 : iReale et al~l l2010h . DEM analyses pro - 



vided ambiguous answers (Schmelz et al 
Several studies ( Feldman et al. 19981 



20091) 



19991: lLandi fc Feldman 1 120081 ) suggest the ubiq- 
uitous presence of isothermal plasma in the quiet 
corona. If confirmed, these results would chal- 
lenge many theoretical models, but the relia- 
bility of the tempera ture diagnostics u sed has 



been questioned (e.g. lLandi et al. 1 120121 ). Part 



of this state of facts is due to technical issues, 
such as t he difficulty to subtract the background 



emission (Terzo fc Realdl2010l: Aschwanden et al 



2008) or the possible spatial and temporal mis- 



Warren 



match of the structures observed in different bands 
and/or different instruments. But the fundamen- 
tal limitations of the DEM inversion clearly play 
a role. 

Despite the many proposed inversion schemes, 
rigorously estimating levels of confidence in the 
various possible solutions given the uncertainties 
remains a major difficulty. In this work we pro- 
pose a strategy to explore the whole parameter 
space, to detect the presence of secondary solu- 
tions, and to compute their respective probabil- 
ities. We illustrate the method by characteriz- 
ing the inversion of simulated observations of the 
AIA telescope. This approach allows us to quan- 
tify the robustness of the inversion by comparing 
the inverted DEM to the input of the AIA simula- 
tions. The main drawback is the limitation to sim- 
ple DEM forms (i.e. described by a small number 
of parameters), but the results provide important 
insights into the properties of more generic DEM 
inversions. The method can also easily be applied 
to any instrument, broad band imagers and spec- 
trometers alike. 

Paper I was dedicated to the analysis of isother- 
mal plasmas, and showed the existence of multi- 
ple solutions in the case of a limited number of 
bands. Notwithstanding, the statistical method 
developed in this work (which provides the respec- 
tive probabilities of each secondary solution), en- 
ables to properly interpret the solutions as consis- 
tent with several plasma temperatures. The use of 
the six AIA coronal bands was shown to increase 
the thermal diagnostic capabilities, providing a ro- 
bust reconstructions of isothermal plasmas. The 
detailed analysis of the squared residuals showed 
that the DEM complexity is limited by the redun- 
dancy of information between the bands. 

We now generalize the method to DEM models 
able to describe a great variety of plasma condi- 
tions, from isothermal to broadly multithermal. A 



2 



summary of the technique is given FigureQ]and its 
implementation is exhaustively described in Paper 
I. The core of our method resides in the proba- 
bilistic interpretation of the DEM solutions. The 
DEM £ 7 solution of the inversion is the one that 
minimizes a criterion C(£) defined as the distance 

observed in iV& bands 
. We write 



between the intensities s 
and the theoretical ones if 1 



« p )-«) 



6=1 

g =ajgminC(£) 
X 2 = minC(£) 



(3) 



where a% is the standard deviation of the un- 
certainties and x 2 represents the residuals. Us- 
ing Monte-Carlo simulations of the systematic and 
random errors, we compute -P(£ 7 |£ p ) the condi- 
tional probability to obtain a DEM £ J knowing the 
plasma DEM £ p . Both types of errors are mod- 
elled as Gaussian random variables with a 25% 
standard deviation (see Section 2.3.3 of Paper I). 
Then, using Bayes' theorem we obtain -P(f p |£ 7 ), 
the probability that the plasma has a DEM £ p 
given the inverted DEM ^ . Using this latter 
quantity it is possible to identify the range or mul- 
tiple ranges of solutions consistent with a set ob- 
servations. Therefore, even if it is not possible 
to alleviate the degeneracy of the solutions, it is 
at least possible to notice and to quantify them. 
But the computation of the probability -P(^ P |C / ) 
is practical only if the space of solutions is lim- 
ited. We therefore restrict the possible solutions 
to simple forms described by a small number a 
parameters: Dirac delta function, top hat, and 
Gaussian. We thus do not to propose a generic 
inversion algorithm, but using this approach, we 
were able to completely characterize the behavior 
of the inversion in well controlled experiments. 

All the results presented in the present paper 
concern multithermal plasmas. In Section [2J we 
analyse the properties of the isothermal solutions 
in order to determine to what extent it is possible 
to discriminate between isothermal and multither- 
mal plasmas. Section [3] is then dedicated to the 
general case of multithermal inversions. Results 
are summarized and discussed in Section 01 



2. Isothermal solutions 

As already mentioned, a recurring question is 
that of the isothermality of solar plasmas . The 



EM loci technique (e.g. Del Zanna et al. I [2002; 



Del Zanna fc Masonll2003l and references therein) 
was originally proposed as an isothermality test. 
The EM loci curves are formed by the set of 
(EM,T e ) pairs for which the isothermal theoret- 
ical intensities exactly match the observations in 
a given band or spectral line. The isothermal hy- 
pothesis can then in principle be ruled out if there 
is no single crossing point of the loci curves. But a 
fundamental ambiguity arises from the inevitable 
presence of measurement errors. 

Under the hypothesis of ideal measurements, 
the observations and the theoretical intensities are 
both reduced to 7°. Hence, the isothermal solu- 
tion for a perfectly isothermal plasma leads to a 
residual \ 2 exactly equal to 0, and the EM loci 
curves all intersect at a common point. But, in 
reality, the presence of errors leads to over or 
under-estimations of the theoretical and observed 
intensities, which conducts to a non-zero \ 2 - I n 
this case, even though the plasma is isothermal, 
the EM loci curves do not intersect at a single 
point, each one being randomly shifted from its 
original position. But even with perfect measure- 
ments, the isothermal hypothesis would also yield 
a residual greater than zero if the observed plasma 
is multithermal, because in this case the I® and 
I^ h are intrinsically different. Therefore, measure- 
ment errors can be incorrectly interpreted as devi- 
ations from isothermality and vice-versa, a variety 
of multithermal plasmas can be statistically con- 
sistent with the isothermal hypothesis. The ques- 
tion is then whether or not a perfectly isothermal 
plasma can be distinguished from a multithermal 
one and if so, under what conditions. 

We thus investigate in this section the isother- 
mal solutions £( s0 obtained from simulated obser- 
vations Ig bs of plasmas of different degrees of mul- 
tithermality. The adopted statistical approach al- 
lows us to quantify to what extent the robustness 
of the inversion is affected as a function of the de- 
gree of multithermality of the plasma. The plasma 
is assumed to have a Gaussian DEM defined in 
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logT e as 

£ p au (EM,T c ,a) =EMN(\ogT e - logT c ), 

Kf ( \ 1 ( x2 \ (4) 

with Al (X) — — = exp — — - 

where EM is the total emission measure, T c is 
the central temperature and a is the width of the 
DEM. It represents a plasma predominantly dis- 
tributed around a central temperature T c and can 
represent both isothermal and very multithermal 
plasmas. 

The theoretical intensities l£ h are limited to the 
case of isothermal solutions, corresponding to a 
DEM 

$ iso (EM, T c ) = EM 5{T e - T c ), (5) 

where <5 is Dirac's delta function and EM is the 
total emission measure. The solution £ J of the 
inversion process then gives the pair (T c 7 , EM 1 ) 
that best explains the observations. The reference 
intensities 7° (see the previous Section), used to 
compute both Ig bs and if 1 have been tabulated 
once and for all for 500 central temperatures from 
logT c (K) = 5 to 7.5, 200 emission measures from 
log.EM(cm -5 ) = 25 to 33, and for the Gaussian 
model 80 DEM widths from to 0.8, expressed in 
unit of logT e (see Section 2.3.2 of Paper I). 

2.1. Ambiguity between uncertainties and 
mult it her mality 

In Figure H the criterion C(EM, T c ) (Equa- 
tion [3]) is plotted for different plasma configura- 
tions. It is represented in gray scal^_] and its ab- 
solute minimum, corresponding to the solution ^ , 
is marked by a white plus sign. The criterion is a 
function of the two parameters EM and T c defin- 
ing the isothermal solutions (Equation ([5])). The 
white curves are the EM loci curves, correspond- 
ing to the minimum of the contribution of each 
band to the total criterion. In all panels, the sim- 
ulated plasmas have DEMs centered on the typical 
coronal temperature T c p = 1.5 x 10 6 K, and a to- 
tal emission measure EM P = 2 x 10 29 cm -5 , an 
active region value guaranteeing signal in all six 
bands (see Section 3 and Figure 3 in Paper I) . The 



X A11 the criteria presented in this paper 
are available in color on line at ftp.ias.u- 
psud.fr/cguennou/DEM_AIA_inversion/. 



top left panel corresponds to an isothermal plasma 
(c = which corresponds to the case presented in 
Section 3.2 of Paper I). Because of systematics, 
Sfc, and instrumental noises, rib, and thus the over 
or under-estimation of the observed and theoret- 
ical intensities I^ bs and 7^ , the EM loci curves 
are shifted up and down along the emission mea- 
sure axis, with respect to the zero-uncertainties 
case (i.e. rib and Sb = 0). Therefore, the curves 
never intersect all six at a single position, thus 
giving x 2 > 0. But multithermality has the same 
effect: the criterion presented on the top right cor- 
responds to a slightly multithermal plasma having 
a Gaussian DEM width of a p = 0.1 logT e . Even 
though in this case uncertainties have not been 
added, the loci curves and the absolute minimum, 
corresponding to the location where the curves are 
the closest together, are also shifted relative to 
each other along the emission measure axis and 
there is no single crossing point. Thus, errors 
and multithermality can both produce compara- 
ble deviations of the observations from the ideal 
isothermal case, hence the fundamental ambiguity 
between the two. The two criteria of the bottom 
panels have been obtained with the same plasma 
Gaussian DEM distribution, but now in presence 
of random perturbations. These two independent 
realizations of the uncertainties illustrate an ex- 
ample of the resulting dispersion of the solutions. 

In Figure El we increased the DEM width to 
a p = 0.3 logTe (top row) and a p = 0.7 logT e 
(bottom row) while keeping the same emission 
measure and central temperature. For each width, 
the left and right panels show two realizations of 
the uncertainties. The corresponding configura- 
tion of the loci curves and thus the value and po- 
sition of the absolute minimum can greatly vary, 
even though the plasma is identical in both cases. 
It also appears that it exists privileged tempera- 
ture intervals where the solutions tend to concen- 
trate. This phenomenon is intrinsically due to the 
shape of the EM loci curves. In the next subsec- 
tion, we characterize the respective probabilities 
of occurrence of the various solutions for a wide 
range of plasma conditions. From the examples 
presented in Figure El we also note that the verti- 
cal spread of the loci curves tends to be larger for 
wider DEMs, which corresponds to larger \ 2 resid- 
uals. In section [2731 we will determine up to what 
DEM width the observations can be considered to 
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be consistent with the isothermal hypothesis. 

2.2. Probability maps 

We consider Gaussian DEM plasmas and scan 
all possible combinations of central temperatures 
and widths used to pre-compute the reference the- 
oretical intensities. The total emission measure is 
kept constant at EM P = 2 x 10 29 cm" 5 . For each 
combination of plasma parameters, 5000 Monte- 
Carlo realizations of the random perturbations rib 
and systematics Sb are obtained, leading to 5000 
estimates of the observed intensities Ib bs {£, P ) m 
each band b. The isothermal least-square inver- 
sion of these 5000 sets of AIA simulated observa- 
tions provides as many solutions £ 7 , allowing an 
estimation of P(£ 7 |£ p ). 

Figures HI [5] and [6] show the resulting temper- 
ature probability maps for a p = 0.1, 0.3 and 
0.7 logT e respectiveljo The isothermal plasma 
case (cr p = 0) is shown in Figure 6 of Paper I. In 
all figures, the probability P{Tl\T p ) is obtained 
by reading vertically panels (a), given the prob- 
ability to find a solution T/ knowing the plasma 
central temperature T p . The probability profiles 
for two specific plasma temperatures, 1.5 x 10 6 and 
7 x 10 6 K, are shown in panels (b) and (c). Using 
Baye's theorem, the probability maps P(T P |T/) of 
panels (e) are obtained by normalizing P(T C 7 |T P ) 
to P(T/), the probability of obtaining T c 7 what- 
ever T P (panels (d)). In the bottom right panels 
(f) and (g), we show two example horizontal pro- 
files giving the probability that the plasma has a 
central temperature T P knowing the inverted tem- 
perature T c 7 . 

In panels (a) and (e) of Figure|4l the most prob- 
able solutions are located around the diagonal. 
However, compared to the case of an isothermal 
plasma (see Figure 6 of Paper I), the distribution is 
wider, irregular and deviations from the diagonal 
greater than its local width are present. As shown 
by panel (d) and the nodosities in the map (a) , the 
unconditional probability to obtain a result T c 7 is 
non uniform, meaning that some inverted temper- 
atures are privileged whereas others are unlikely. 
Comparing with Figure 6 of Paper I, profile (b) 
shows that the probability of secondary solutions 



2 The probability maps for 80 widths from to 
0.8 logTe are available in color on line at ftp.ias.u- 
psud.fr/cguennou/DEM_AIA_inversion/. 



at T C P = 1.5 x 10 6 K is increased with respect to 
the isothermal case. The apparition of these two 
solutions are illustrated in the bottom row of Fig- 
ure [5J The bottom right panel corresponds to a 
realization of the errors yielding a solution close 
to the diagonal, while the bottom left panel of the 
same Figure illustrates a case where the absolute 
minimum of the criterion is located at low temper- 
ature. Using the map of P(T P \T*) it is however 
possible to correctly interpret the low tempera- 
ture solutions as also compatible with 1.5 x 10 6 K 
plasma (profile (g)). 

In Figure 03 the plasma DEM width is in- 
creased to a p = 0.3 logT e . As a result, the 
above-described perturbations with respect to the 
isothermal plasma case are amplified. The diag- 
onal structure has almost disappeared, with dis- 
continuities and reinforced and more diffuse no- 
dosities. Multiple solutions of comparable prob- 
abilities are present over large ranges of plasma 
temperatures and consequently, the estimated T 7 
can be very different from the input T p . For 
example, panel (c) shows that for a 7 x 10 6 K 
plasma, the most probable T c 7 is either 1.6 x 10 5 
or 3 x 10 5 K. The unconditional probability P(P C 7 ) 
of panel (d) is very non uniform, some ranges 
of estimated temperatures being totally unlikely 
(e.g. Tl = [1.5 x 10 6 ,4 x 10 6 ] K) while oth- 
ers are probable for large intervals of T p (e.g. 
Tl = 3 x 10 5 K or 10 6 K). However, despite the 
jaggedness of P(T C 7 |T C P ), the map of P(T C P |T C 7 ) can 
once again help properly interpret the result of the 
inversion. For example, profile (g) shows that for 
Tl = 1.5 x 10 6 K, the distribution of T p is dis- 
tributed around T c p = 10 7 K, which is exactly the 
plasma temperature that can yield an inversion 
at T/ = 1.5 x 10 6 K (see panel (a)). Panel (f), 
providing the probability distribution T P know- 
ing that the inversion results is T 7 = 7 x 10 6 K, 
exhibits a broad probability distribution around 
T p = 1.5 x 10 7 , showing that the plasma temper- 
ature thus deduced is very uncertain. This is to 
be compared to the 0.05 logT c p temperature res- 
olution of the isothermal case (see section 3.2 of 
Paper I). 

As the DEM becomes even larger, the im- 
pact on the robustness of the inversion becomes 
greater. At a p = 0.7 logP e , the probability map 
P(P C 7 |T P ) of Figure [6] (a) and the correspond- 
ing probability P(T 7 ) show clearly two privileged 
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solutions: T/ = 10 6 and 3 x 10 5 K. The es- 
timated isothermal temperatures are always the 
same for any T p , as illustrated by panels (b) and 
(c). Therefore, the inversion results contains 
no information on the plasma central temperature 
T c p . This is illustrated by the lack of structure in 
the probability map P{T c p \T,f) of panel (e). Pro- 
file (g) shows that for T/ = 1.5 x 10 6 K, the dis- 
tribution of T p extends all over the temperature 
range. 

2.3. Interpretation 

We have shown that as the width of the plasma 
DEM increases, multiple solutions to the isother- 
mal inversion appear. This phenomenon has been 
alread y mentioned by iPatsourakos fc KlimchukT 



(|2007l ) using triple-filter TRACE data. After a 
proper treatment of the uncertainties, the authors 
found that their observations of coronal loops were 
consistent with both a high (w 1.5 x 10 6 K) and 
a low (w 5 x 10 5 K) isothermal plasma tempera- 
ture. They correctly concluded that without a pri- 
ori knowledge of the physical conditions in these 
loops, they could not rule out the cool plasma so- 
lutions. Even though we used six bands, multiple 
solutions appear anyway with increasing plasma 
width. In addition, as we have seen in Paper I, 
multiple solutions can exist even with an isother- 
mal plasma if only a limited number of bands is 
available. This is another illustration of the simi- 
lar effects of errors and multithermality. 

The isothermal temperature solutions become 
progressively decorrelated from the plasma cen- 
tral temperature as the width of the DEM in- 
creases. For very large DEMs (Figure [B]) , the in- 
version process yields exclusively either 3 x 10 5 K 
or 10 6 K whatever the plasma T P . These two 
temperatures correspond to the preferential loca- 
tions of the minima shown in the criteria of Fig- 
ure [3] This i s a generalizat i on of t he phenomenon 
analyzed by IWeber et al.l (l2005h in the simpler 
case of the TRACE 19.5 over 17.3 nm filter ra- 
tio. The authors showed that in the limit of an 
infinitely broad DEM, the band ratio tends to a 
unique value equal to the ratio of the integrals of 
the temperature response functions. Furthermore, 
they showed that as the width of the DEM in- 
creases, the temperature obtained from the band 
ratio becomes decorrelated from the DEM central 
temperature. We have found a similar behavior in 



the more complex situation of six bands. This is 
not however an intrinsic limitations of AIA. We 
can predict that the same will occur with any 
number of bands or spectral lines. Indeed, for an 
infinitely broad DEM, since the observed inten- 
sities are equal to the product of the total emis- 
sion measure by the integral of the response func- 
tions (Equation [2|), they are independent from the 
plasma temperature. Therefore, the inversion will 
yield identical results for any plasma temperature 
T c p , whatever the number of bands or spectral 
lines. 

2.3.1. Defining isothermality 

As already noted in section [2TT1 and in Figure [2 
to larger DEM widths correspond larger squared 
residuals. From Paper I the distribution of resid- 
uals to be expected for an isothermal plasma is 
known. Examining then the residuals for the so- 
lutions given in the probability maps of section [2~2l 
all the solutions may not all be statistically con- 
sistent with the isothermal hypothesis. We will 
thus analyze the distribution of residuals to define 
rigorously a test of the adequacy of the isothermal 
model used to interpret the data. 

Because both the random and systematic er- 
rors ni, and Sb have been modeled by a Gaussian 
random variable, if the DEM model used to in- 
terpret the data can represent the plasma DEM, 
the residuals are equal to the sum of the square 
of six normal random variables (see Equation [3]) . 
Since we adjust the two parameters EM and T c , 
the residuals should thus behave as a degree four 
X 2 distribution. 

Figure [7] shows the distribution of the squared 
residuals for all the 80 DEM widths considered in 
the simulations. The shades of grey in the top 
panel correspond to the probability to obtain a 
given x 2 value (abscissa) as a function of the width 
of the plasma DEM (ordinate). In the bottom 
panel, four profiles give the distribution of squared 
residuals obtained for an isothermal plasma (thin 
dotted line) and for the three DEM widths dis- 
cussed in section 12.21 a p = 0.1 (thick dotted 
line), 0.3 (thick dashed curve) and 0.7 logT e (thick 
dash-dotted curve), corresponding to the white 
horizontal lines on the top panel. The theoretical 
X 2 distributions of degree three (thin solid curve) 
and four (dashed curve) are also plotted. 
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If the plasma is isothermal, the distribution of 
the residuals is slightly shifted toward a degree 3 
X 2 instead of the expected degree four. In paper I, 
we interpreted this as a correlation between the six 
AIA coronal bands. The distributions of residuals 
progressively depart from the isothermal case as 
the DEM width of the simulated plasma increases. 
The distributions become broader and their peaks 
are shifted toward higher values, forming the diag- 
onal structure in the top panel of Figure [JJ This 
behavior stops around a p — 0.4 logT e . Above, 
the peaks of the distributions are shifted back to- 
ward smaller values and remains constant. As the 
DEM becomes wider, the simulated observations 
become independent from the plasma parameters, 
and all inversions tend to give the same solution 
and the same residuals. 

These distributions of residuals provide a ref- 
erence against which to test the pertinence of the 
isothermal model. The isothermal hypothesis can 
for example be invalidated for the solutions bi- 
ased towards T* » 1 MK and corresponding to 
very multithermal plasmas (e.g. a p = 0.7, see 
section 12.21 and Figure |6]). Indeed, for a given 
inversion and its corresponding residual, the top 
panel of Figure [7J gives the most probable width 
of the plasma, assuming it has a Gaussian DEM. 
Let us assume that an isothermal inversion re- 
turns a residual equal to 5. Analyzing the his- 
togram corresponding to the bottom row of the 
top panel of Figure [JJ (er p =0), we can show that 
an isothermal plasma has 68.2% chance to yield 
X 2 < 5. This residual can therefore be consid- 
ered consistent with an isothermal plasma. But 
reading the plot vertically, we see that the proba- 
bility P(x 2 = 5) is greater for multithermal plas- 
mas and peaks around a p = 0.12, which is thus 
in this case the most probable Gaussian width. 
For larger residuals the situation is more complex 
because several plasma widths can have equally 
high probabilities. Past \ 2 = 3.5 the plasma has 
a higher probability to be Gaussian than isother- 
mal and a Gaussian inversion is required to prop- 
erly determine its most probable width and central 
temperature. 

Figure I7J is a generalization of th e results ob- 
tained bv Eandi fc Klimchuk I (|2010f ). It is equiv- 



alent to their Figure 2, identifying our x 2 to their 
-Fmin and their solid and dashed lines to respec- 
tively the peak and the contours at half maximum 



of our x 2 distribution as a function of a p . Their 
Figure 2 was computed using 13 individual spec- 
tral lines for a 1 MK plasma and extends only up 
to a p rj 0.2, while our Figure [7J was computed 
for the 6 AIA bands over a wide range of central 
temperatures and for widths up to a p = 0.8. De- 
spite these differences, the two figures exhibit the 
same global behavior. Indeed, for any number of 
spectral lines or bands, the residuals of the isother- 
mal inversion tend to increase as the width of the 
plasma DEM increases. 

3. Multithermal solutions 

We now focus on multithermal solutions. The 
ability of AIA to reconstruct the DEM given the 
uncertainties is evaluated, and the probability 
maps associated to all parameters are computed, 
allowing to take into account all the Gaussian 
DEMs consistent with the simulated observations. 
Both cases of consistent and inconsistent DEMs 
models between the simulations and the inversion 
assumptions are examined. After the previous sec- 
tion on isothermal solutions, this generalizes the 
study of the impact of a wrong assumption on the 
DEM shape. 

As in the previous section, the simulated ob- 
servations remain Gaussian, but we now consider 
Gaussian solutions, i.e. theoretical intensities I^ 1 
tabulated for the Gaussian DEM model. The 
model can thus in principle perfectly represent the 
plasma conditions. 

3.1. Three-dimensional criterion 

Investigating now multithermal Gaussian solu- 
tions, the least-square criterion given by Equa- 
tion ([3]) has thus three dimensions C(£ gau ) = 
C(EM,T c ,a). This three dimensional parameter 
space is systematically scanned to locate the the- 
oretical intensities describing the best the simu- 
lated observations. Figure [5] shows this criterion 
for three cases illustrating different degrees of mul- 
tithermality. In each row, from top to bottom, the 
simulated plasma has a DEM width a p of 0.1, 0.3 
and 0.7 logT e respectively, centered on the tem- 
perature T c = 1.5 x 10 6 K. On the left panels, the 
background image represents C(EM,T c ,a I ), the 
cut across the criterion in the plane perpendicu- 
lar to the DEM width axis at the width a 1 cor- 
responding to the absolute minimum (white plus 
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sign). The curves are the equivalent of the loci 
emission measure curves in a multithermal regime: 
for each band b and for a given DEM width er p , 
they represent the loci of the pairs (EM,T C ) for 
which the theoretical intensities Ijj h are equal to 
the observations I£ bs . As a varies, they thus de- 
scribe a loci surface in the three-dimensional cri- 
terion. The difference with the isothermal loci 
curves is that the theoretical intensities ijj' 1 have 
been computed for the multithermal case (i.e. con- 
sidering a Gaussian DEM), and thus the param- 
eter a must be now considered. The right pan- 
els of Figure E] display C{EM I ,T c ,a), the cuts 
across the criterion in the plane perpendicular to 
the emission measure axis at the EM 1 correspond- 
ing to the absolute minimum of the criterion (also 
represented by a white plus sign). 

The impact of multithermality on the criterion 
topology is clearly visible on the loci curves, in- 
ducing a smoothing of the curves as the multi- 
thermality degree increases. Indeed, the intensi- 
ties I® as a function of the central temperature 
can be expressed as the convolution product be- 
tween the instrument response functions and the 
DEM (see Section 2.3.2 of Paper I). Therefore, 
the reference intensities I®(EM,T C , a) computed 
for multithermal plasmas are equal to the isother- 
mal ones smoothed along the T c axis. As a result, 
the criterion exhibits a smoother topology and the 
minimum areas become broader and smoother as 
the multithermality degree increases, introducing 
more indetermination in the location of the ab- 
solute minimum. For each row, the realization 
of the uncertainties rib and Sb yields a solution 
that is close to the simulation input. The cuts are 
therefore made at similar locations in the criterion 
in order to best illustrate the modification of its 
topology. We will now analyse to what extent the 
distortion of the criterion affects the robustness of 
the inversion. 

3.2. Probability maps 

In the case of DEM models defined by three pa- 
rameters, and since the plasma EM P is fixed, the 
probability matrices P{EM 1 T c 7 , a r \EM F ', T P , a p ) 
resulting from the Monte-Carlo simulations have 
five dimensions. In order to illustrate the main 
properties of these large matrices, we rely on com- 
binations of fixed parameter values and summa- 
tion over axes. 



The associated probability maps are displayed 
in Figure |9] for simulated plasmas characterized 
by g p = 0.1 (top panel) and a p = 0.7 logT e 
(bottom paneljj. The probability maps are repre- 
sented whatever the emission measure and DEM 
width obtained by inversion, i.e. the probabili- 
ties P(EM I ,T I ,a I \EM p ,T p ,a p ) are integrated 
over EM 1 and a 1 . In case of a narrow DEM dis- 
tribution (er p = 0.1, top panels) the solutions are 
mainly distributed along the diagonal, with some 
secondary solutions at low probabilities. The ac- 
curacy of the determination of the central tem- 
perature T p is improved compared to the isother- 
mal inversion of the same plasma as shown in Fig- 
ure |H the diagonal is more regular and P(T/) is 
more uniform. However, increasing the width of 
the DEM of the simulated plasma reduces the ro- 
bustness of the inversion process. In the bottom 
panels (a p = 0.7), we observe a distortion and an 
important spread of the diagonal. In particular, 
the horizontal structure in panel (a) and the con- 
sequent peak of P(T/ ) tell that the estimated tem- 
perature is biased toward T 1 ~ 1 MK, especially 
for plasma temperatures T p < 2 MK. This is to 
be compared with the privileged isothermal solu- 
tions in Figure O and the same reasoning applies. 
For very broad DEMs, because of the smoothing of 
the If}(EM,T c ,a) along the T c axis, the observed 
intensities are only weakly dependent on the cen- 
tral temperature, and thus all inversions tend to 
yield identical results. As a consequence, the prob- 
ability map P{T P \T') shows that over the whole 
temperature range considered there is a large un- 
certainty in the determination of T P . The smooth- 
ing of the 1° (EM, T c , a) is due to the width of the 
plasma DEM and not to the properties of the re- 
sponse functions. It is therefore to be expected 
that the uncertainty in the determination of T p 
persists even if using individual spectral lines. 

In the same way, the conditional probabilities 
P(o-V P ) and P(cr P |cr / ) of the DEM width are 
displayed on Figure [TU] The probabilities are 
represented whatever the estimated temperature 
Tl and emission measure EM 1 , for a simulated 
plasma having a DEM centered on T P = 1 MK. 
The conditional probability P(<7 7 |<7 P ) on panel (a) 
exhibits a diagonal that becomes very wide for 



3 The probability maps for 80 widths from to 
0.8 logT e are available in color on line at ftp.ias.u- 
psud.fr /cguennou/DEM_AI Ajnversion / . 
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large a p and from which another broad branch 
bifurcates at a 1 = 0.15 logT e . This branch is the 
analog of those observed on the temperature axis 
(top panels of Figure ^ . For plasmas having a 
DEM width greater than a p = 0.3 logT e , the esti- 
mated width a 1 is decorrelated from the input a p . 
Profiles (b) and (c) provide the conditional prob- 
abilities of a 1 for plasma DEM widths a p = 0.2 
and 0.6 log T e respectively. As shown by panel (d), 
the unconditional probability P(a ! ) to obtain a 1 
whatever T c 7 for this 10 6 K plasma is biased to- 
wards a 1 = 0.12 logT e . Exploring the full prob- 
ability matrix P{EM I , T/, cr 7 |£M p , T c p , a p ), we 
discovered that if the plasma DEM is broad, the 
small width solutions (« 0.12 logT e ) are the ones 
that were also biased towards a central temper- 
ature of 10 6 K in Figure El This is caused by 
the shape of the temperature response functions 
Rb(T e ). A minimum is formed in the criterion 
C(T c ,a) around (T c = 1 MK, a = 0.12 logT e ) 
that tends to be deeper than the other local min- 
ima. Since for a very broad DEM plasma the sim- 
ulated observed intensities I^ hs are almost always 
similar, most of the random realizations conduct 
to this deeper minimum and to the formation of 
the narrow solutions branch in P(a I \a p ) panel 
(a). Therefore, very multithermal plasmas tend 
to be systematically measured as near-isothermal 
and centered on 10 6 K. 

Reading horizontally in panel (e) the inverse 
conditional probability P(a p \a I ), a large range of 
DEM widths a p is consistent with the estimation 
a 1 . Both plots (f) and (g), representing two cuts 
at = 0.2 and 0.6 logT c , can be used to deduce 
the most probable plasma DEM with a p . How- 
ever, at a 1 — 0.6 logT e , almost all plasma widths 
have significant probabilities, restricting consid- 
erably the possibility of inferring relevant DEM 
properties. The situation improves for smaller 
widths, as shown in panel (g), even though the 
probability distribution is still broad. For narrow 
DEM distributions (er p < 0.2), the width of the 
distribution decreases to about 0.15 logT e , which 
thus represents the width resolution limit. 

The analysis of the probability maps demon- 
strates that the robustness of the inversion is sub- 
stantially affected by the degree of multithermality 
of the observed plasma. Furthermore, as we al- 
ready noted, the simulations presented have been 
made in a favorable configuration where signifi- 



cant signal is present in all six bands. For narrow 
plasma DEMs (a p < 0.15 logT e ), the six AIA 
coronal bands enable an unambiguous reconstruc- 
tion of the DEM parameters within the uncertain- 
ties. The precision of the temperature and DEM 
width reconstruction is then given by the widths 
of the diagonals in Figures |9] and [TQl For the tem- 
perature, that width is consistent with the reso - 
lution of [0.1,0.2] logT e given bv Ijudge I (|2010l ). 
However, both the accuracy and the precision of 
the inversion decrease as the multithermality de- 
gree of the simulated plasma increases, a wider 
range of solutions becoming consistent with the 
observations. In the case of a very a large DEM 
plasma, the solutions are skewed toward narrow 
1 MK Gaussians. The isothermal solutions bi- 
ased towards 1 MK (section 12 . 3[) could be inval- 
idated on the basis of their correspondingly high 
residuals (section [2.3. 1|) . But the biased Gaussian 
solutions are by definition fully consistent with 
a Gaussian plasma. T his result generalizes that 
of Weber et al. I (2JD05), and since our method en- 



sures that the absolute minimum of the criterion 
is found, it is a fundamental limitation and not an 
artefact of the minimization scheme. Instead of 
being evidence for underlying physical processes, 
the recurrence of common plasma properties de- 
rived from DEM analyses may be due to biases in 
the inversion processes. However, using the type 
of statistical analysis presented here, it is possible 
to identify theses biases and correct for them. 

3.3. Residuals and model testing 

The distribution of the sum of squared residuals 
is displayed in the top panel of Figure [TT] as a func- 
tion of the plasma DEM width, a p , of the simu- 
lated plasma. The three white horizontal lines rep- 
resent the locations of the cuts at a p = 0.1 (dot- 
ted line), 0.3 (dashed line) and 0.7 logT e (dotted 
dashed line) shown in the bottom panel. Unlike 
for the isothermal solutions (Figure [3]) , the distri- 
butions are similar for all plasma widths and re- 
semble a degree 3 \ 2 distribution (thin solid line) . 
The most probable value is about 1 and 95% of the 
residuals are between and 10. Any DEM inver- 
sion yielding a x 2 smaller than 10 can thus be con- 
sidered consistent with the working hypothesis of 
a Gaussian DEM plasma. It does not mean how- 
ever that a Gaussian is the only possible model, 
but that it is a model consistent with the observa- 
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tions. Since we perform a least square fit of the six 
values Iobs_jth by three p aram eters (EM, T c , a), 
this behavior is expected. The small correlation 
between the residuals observed in the isothermal 
case disappears with increasing DEM width. This 
can be explained by the smoothing of the criterion 
(see Figured]) that reduces the directionality of its 
minima. 

For model testing, a reduced chi-squared x 2 e d = 
X 2 /n, where n is the number of degree of freedom, 
is sometimes preferred over a regular \ 2 test, for 
it has the advantage of being normalized to the 
model complexity. Usually, n is the number of 
observations minus the number of fitted param- 
eters. This assumes that the measurements are 
independent which is what we wanted to test. In 
practice, a small correlation between the residuals 
is found in the isothermal case, but the effect is 
small and a reduced % 2 can be used. It should 
be noted that the residuals have a non negligible 
probability to be greater than the peak of the cor- 
responding \ 2 distribution. For plasmas having 
Gaussian DEMs for example, the residuals have 
about a 43% chance to be greater than 3 (Fig- 
ure [XT]) . This implies that seemingly large resid- 
uals are not necessarily a proof of the inadequacy 
of the chosen DEM model. The working hypothe- 
sis can only be invalidated if it can be shown that 
another model has a greater probability to explain 
the obtained residuals. This was for example the 
case in Section 12.3.11 where we have shown that 
if an isothermal inversion gives a squared residual 
greater than 3, the plasma DEM is more proba- 
bly Gaussian than isothermal. We now explore a 
similar situation for Gaussian inversions. 

Indeed, in reality the DEM shape is not known 
and it is interesting to test if a Gaussian is a per- 
tinent model or if AIA has the capability to dis- 
criminate between different models. For this pur- 
pose, still considering Gaussian DEM solutions, 
the simulations of the observations are now per- 
formed using a top hat DEM distribution defined 
as 

& ot (T e ) =EM n Tc , CT (T e ), 

'i if|log(T e )-log(T c )|<f 



with n TcjCT (T e ) = i* 

else 



(6) 



Like a Gaussian, this parametrization can repre- 



sent narrow and wide thermal structures. 

Figure [12] gives the associated temperature 
probability maps. The top panels correspond to 
simulated plasmas with a top hat distribution of 
width <r p = 0.1. Most of the solutions are con- 
centrated around the diagonal, even though the 
robustness is somewhat affected for temperatures 
in the range 5 x 10 5 < T c p < 10 6 K, where low 
probability secondary solutions exist. These two 
plots are very similar to the top panels of Figure [5] 
that give the probability of the Gaussian solutions 
to a Gaussian plasma. Therefore, even though the 
DEM assumed for the inversion is different from 
that of the plasma, the central temperature of the 
top-hat is determined unambiguously. 

The picture is different for wide plasma DEM 
distributions. In the bottom panels of Figure [121 
with a p = 0.7, the diagonal has become very wide, 
structured and does not cross the origin any more. 
Some plateaus appears, meaning that for some 
ranges of plasma temperature T c p , the tempera- 
ture T/ provided by the Gaussian solutions will be 
invariably the same. For 5 x 10 5 < T c p < 3 x 10 6 K 
for example, a constant solution T J = 1 MK ap- 
pears, as shown also by P{T^) (bottom panel (b)). 
As a result, the inverse conditional probability 
map P(T P |T C 7 ) (bottom panel (c)) indicates that 
for an inversion output of 1 MK there is a large 
indetermination on the central temperature. The 
behavior of the solutions is thus globally equiva- 
lent to that described in Section 13.21 for the inver- 
sion of Gaussian DEM plasmas. 

The observed distribution of the sum of the 
squared residual is very similar to those obtained 
with consistent Gaussian DEM models and is close 
to 3 x 2 ( see Figure [5]). If the working hypothe- 
sis were a top hat DEM, thus consistent with the 
plasma, the residuals would also be close to a de- 
gree 3 x 2 : as f° r an Y DEM model described by 
three parameters (see discussion above). The so- 
lar plasma has no reason to really have a top hat 
DEM. However, this numerical experiments shows 
that, using AIA data only and the x 2 test, the 
discrimination between two very different multi- 
thermal DEMs is practically impossible. 

4. Summary and discussion 

In this work we described a complete charac- 
terization of the statistical properties of the DEM 
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inversion, rigorously treating both systematic and 
random errors. The developed methodology has 
been illustrated in the specific case of the AIA 
telescope, but the technique is generic and can 
be applied to any other instrument, spectrome- 
ters as well as imaging telescope. By restricting 
ourselves to parametric DEMs, we could analyze 
in details what occurs during the inversion pro- 
cess, and could therefore point out the fundamen- 
tal difficulties involved in the DEM reconstruc- 
tion. Even though only a few simple DEM distri- 
butions have been studied, important and generic 
conclusions regarding the robustness of the inver- 
sion problem have been reached. The method can 
be applied to other forms of DEMs as long as they 
can be defined by a small number of parameters. 

Our technique provides new tools to facilitate 
the interpretation of the DEM inversion. By com- 
puting the -P(£ 7 |£ P ) probability matrices the ro- 
bustness can be evaluated, secondary solutions 
can be detected and their probabilities quantified. 
Since we do not know whether the systematic er- 
rors are over or under-estimated, their random- 
ization in the computation of the inverse condi- 
tional probability P(£ p |£' r ) ensures that all possi- 
ble DEMs consistent with the solution £ J given the 
uncertainties are taken into account. The specific 
location of secondary solutions in the space of the 
possible DEMs £ can be interpreted by analyzing 
the topography of the criterion C(£), which is in- 
trinsically linked to the shapes of the temperature 
response functions. The simulated distributions of 
the sum of the squared residuals provide references 
against which to asses the pertinence of the DEM 
model used to perform the inversion. Applying 
the technique on real data, and applying the tools 
presented here on the results, the interpretation of 
DEM inversions can be facilitated. 

By investigating the isothermal solutions, we 
have shown that the six AIA coronal bands greatly 
enhance the robustness of the DEM reconstruction 
compared to previous EUV imagers, even though 
the temperature resolution varies between 0.03 
and 0.11 logT c . This resolution depends on the 
level of uncertainties assumed (see Section 3.2 in 
Paper I). Because of the ambiguity between mea- 
surement errors and deviation from isothermality, 
a variety of plasma conditions is consistent with 
the isothermal hypothesis. As the degree of mul- 
tithcrmality of the observed plasma increases, the 



isothermal solutions become progressively decorre- 
lated from the input. We have demonstrated that, 
in the limit of very broad multithermal plasmas, 
the isothermal solutions are biased toward two 
preferred solutions: 3 x 10 5 or 10 6 K depending on 
the systematics. But the analysis of the behavior 
of the residuals allows to measure the adequacy of 
the isothermal hypothesis. Past a critical value of 
the residuals, a multithermal DEM is more proba- 
ble than an isothermal plasma. That critical value 
has to be estimated from Monte-Carlo simulations 
for the considered set of bands or spectral lines. 
Our procedure provides a rigorous quantification 
of the classical EM loci isothermal test. 

Investigation of the multithermal Gaussian so- 
lutions shows that AIA has the ability to recon- 
struct simple DEM forms. However, the greater 
the width of the plasma DEM, the lower the ac- 
curacy and precision on the determination of the 
DEM parameters. We argued that this would be 
the case also if using individual spectral lines in- 
stead of narrow bands. In the extreme case of very 
broad DEMs distributions the quality of the inver- 
sion is strongly affected, because of the combined 
effects of a smoothed criterion and simulated ob- 
servations weakly dependent on the plasma pa- 
rameters. The AIA Gaussian solutions are bi- 
asS.P. acknowledges the support from the Belgian 
Federal Science Policy Office through the inter- 
national cooperation programmes and the ESA- 
PRODEX programme and the support of the In- 
stitut d'Astrophysique Spatiale (IAS). F.A. ac- 
knowledges the support of the Royal Observatory 
of Belgium. The authors would like to thank J. 
Klimchuk for fruitful discussions and comments. ed 
towards a near isothermal plasma (a 1 = 0.15) cen- 
tered on T c — 1 MK. Furthermore, the analysis of 
the distribution of the residuals has revealed that 
it is difficult to discriminate between different mul- 
tithermal DEMs models using AIA alone. 

These results can have far reaching implications 
for the determination of the thermal structure of 
the solar corona. Indeed, the thermal structure of 
coronal loops for example is still under debate, dif- 
ferent proposed heating scenarios leading to oppo- 
site DEM properties. For example, the nanoflares 
theor y fe.g. |Parkerlll983tlCargill1ll994tlK"limchuk 



2006), suggesting that the stored energy is then 
released impulsively over small scales, predicts a 
wide thermal structure at a given time, in case of 
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sparse events inside a single strand. Conversely, 
the steady heating processes should produce loops 
with narrow DEMs. Our analysis revealed that 
even if the residuals support the pertinence of 
the solution, the DEM analysis of AIA observa- 
tions can lead to the erroneous conclusion that 
the plasma is nearly isothermal while it is in fact 
broadly multithcrmal. Different biases can be ex- 
pected for different sets of bands or spectral lines, 
different instruments or atomic physics data. 

The methodology developed in this series of pa- 
pers can also be used to establish the optimal set 
of bands needed to maximize the robustness of the 
DEM inversion, for broad band instrument as well 
as spectrometers. Performing this type of prepara- 
tory simulations can thus be very useful to opti- 
mize the design of future solar instruments if they 
are to provide reliable temperature diagnostics. 
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national cooperation programmes and the ESA- 
PRODEX programme and the support of the In- 
stitut d'Astrophysiquc Spatialc (IAS). F.A. ac- 
knowledges the support of the Royal Observatory 
of Belgium. The authors would like to thank J. 
Klimchuk for fruitful discussions and comments. 
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Fig. 1. — Summary of the technique used in this work. The AIA temperature response functions are 
computed using CHIANTI 7.0. The simulated observations Ig bs are calculated assuming a simple DEM model 
£ p (isothermal, Gaussian or top hat), and adding the instrumental random perturbations and systematic 
errors. The theoretical intensities ij h are computed assuming a DEM model £, for a large range of parameters. 
The inverted DEM £ 7 is evaluated by least-square minimization of the distance between the theoretical 
intensities and the simulated observations. Using a Monte-Carlo scheme, TV realizations of the random and 
systematic errors rib and Sb are computed. For a given set of parameters (EM P , T c p , <r p ) of the plasma DEM 
£ p , the quantity P(£ 7 |£ p ) is then evaluated. The probability P(£ P |£ 7 ) is then derived using Bayes' theorem. 
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Fig. 2. — Equivalence between uncertainties and multi-thermality. The criterion C(EM,T C ) given by Equa- 
tion [3l represented in gray scale and its absolute minimum, marked by a white plus sign. The EM loci curves 
for each of the six AIA coronal bands are also represented (white lines) , superimposed to the criterion. The 
location of the absolute minimum provides the isothermal solution £ for a plasma having DEMs £ p centered 
on T c p = 1.5 MK, and a total emission measure of 2 x 10 cm . Top left: isothermal plasma having a DEM 
£ p = £ lso in presence of random and systematic errors. Top right: Multithermal plasma having a Gaussian 
DEM £ p = £3 auss with a width a p = 0.1 logT e , without uncertainties. Uncertainties and multithermality 
both produce comparable deviations of the EM loci curves. Bottom: Same as top right, but for two different 
realizations of the random and systematic errors rib and Sb- 
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Fig. 3. — Isothermal solutions for different thermal widths. Same as Figure [2] for larger DEMs plasma 
widths. In the top row the plasma has a Gaussian DEM a p = 0.3 logT e , and a p — 0.7 logT e in the bottom 
row. For each thermal width, two different realizations of the random and systematic errors are presented, 
showing that the location of the absolute minimum can greatly vary. 
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P(T/\ 1.5MK) P(Tj\ 7MK) 



Fig. 4. — Maps of probabilities considering a Gaussian DEM plasma £ p = £ 9tm having a narrow thermal 
distribution of a p = 0.1 log T e , obtained by 5000 Monte-Carlo realizations of the random and systematics 
errors rn> and s&, and investigating the isothermal solutions, (a): Probability map P(T*\T P ), vertically 
reading. The central temperature resulting from the inversion is presented whatever the total emission 
measure EM 1 . (b) and (c): Probability profiles of T 1 for plasma central temperatures T c p = 1.5 x 10 6 and 
7 x 10 6 K (vertical lines in panel(a)). (d): Total probability to obtain T/ whatever T p (see Section Q] and 
Section 2.1 of paper I), (e): Vice-versa, probability map P(T P \T£), horizontally reading, inferred by means 
of Bayes' theorem using 

p { T c\ T c) and P ( T c)- (I) and (9) : Probability profiles of T P knowing that the 
inversion result is, from top to bottom, 7 x 10 6 and 1.5 x 10 6 K. 
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P(T/\ 1.5MK) P(Tj\ 7MK) 



Fig. 5. — Same as Figure 0] but with a plasma DEM width is increased to a p = 0.3 logT e . Many pertur- 
bations, already visible in the Figure 0] are amplified: the distribution is wider, irregular and the diagonal 
structure disappeared (see panel (a)). The presence of multiples solutions of comparable probabilities is 
increased for a large range of plasma temperatures T^, leading to very different estimated T/ from the 
input T c p . The probability map P(T c 7 |T ( f ) can help us to properly interpret the inversion result, taking into 
account the secondary solutions and providing their respective probability. 
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P(Tj\ 1.5MK) P(Tj\ 7MK) 



Fig. 6. — Same as Figure 21 but in case of the plasma DEM width is increased to a p = 0.7 logT e . The 
impact on the robustness of the inversion is clearly increased in this case, showing two privileged isothermal 
solutions T/, totally decorrelated from the input . As a result, no information regarding the central 
temperature of the DEM can be extracted from the inversion. 
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Fig. 7. — In the top panel, the distributions of the sum of the squared residuals corresponding to isothermal 
inversions of Gaussian DEM plasmas, as a function of the DEM width a p . Higher probabilities correspond 
to lighter shades of grey. The bottom panel shows cuts at a p = (isothermal), 0.1, 0.3 and 0.7 logT e along 
with theoretical \ 2 distributions of degree 3 and 4. 
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Fig. 8. — Criterion in case of Gaussian multithermal DEM inversions. The three-dimensional criterion is 
a function of EM P , T c p , and a p . The simulated plasma has a Gaussian DEM centered on 1.5 MK and a 
different width for each row. From top to bottom: a p = 0.1, 0.3 and 0.7 logT e . The superimposed curves 
on the left panels represent the equivalent of the EM loci emission measure curves in a multithermal regime 
(see the detailed description in Section l3~Tj) . 
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Fig. 9. — Same as figure |4] but investigating now the multithermal solutions. The simulated observations 
are computed for DEM widths of <r p = 0.1 (top row) and 0.7 logT e (bottom row). The probabilities 
are represented here whatever the emission measure EM 1 and the Gaussian width a 1 returned by the 
minimization scheme. In the case of low degree of multithermality plasmas, the solutions remain distributed 
around the diagonal. However, for broad DEM distributions, the solutions appear to be biased toward 
~1 MK, especially for plasmas exhibiting central temperature lower than 2 MK. 
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Fig. 10. — Maps of probabilities for the DEM width for a simulated plasma having a Gaussian DEM 
centered on T c p = 1 MK. Probabilities are represented whatever the emission measure EM 1 and the central 
temperature T/ obtained. The solutions appear to be biased toward a 1 ~ 0.12 log T e whatever the plasma 
widths, even if the situation improves for smaller widths. 
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